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Abstract 

We use Metropolis Monte Carlo and umbrella sampling to calculate the free energies of inter- 
action of two methane molecules and their charged derivatives in cylindrical water-filled pores. 
Confinement strongly alters the interactions between the nonpolar solutes, and completely elimi- 
nates the solvent separated minimum (SSM) that is seen in bulk water. The free energy profiles 
show that the methane molecules are either in contact or at separations corresponding to the diam- 
eter and the length of the cylindrical pore. Analytic calculations that estimate the entropy of the 
solutes, which are solvated at the pore surface, qualitatively explain the shape of the free energy 
profiles. Adding charges of opposite sign and magnitude 0.4e or e (where e is the electronic charge) 
to the methane molecules decreases their tendency for surface solvation and restores the SSM. We 
show that confinement induced ion-pair formation occurs whenever Ib/D ~ 0(1), where Ib is the 
Bjerrum length, and D is the pore diameter. The extent of stabilization of the SSM increases 
with ion charge density as long as Ib/D < 1. In pores with D < 1.2 nm, in which the water is 
strongly layered, increasing the charge magnitude from 0.4e to e reduces the stability of the SSM. 
As a result, ion-pair formation which occurs with negligible probability in the bulk, is promoted. 
In larger diameter pores that can accomodate a complete hydration layer around the solutes, the 
stability of the SSM is enhanced. 
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Introduction: 



Confinement effects on the properties of water and water-mediated interactions between 
solutes are receiving increasing attention because of their relevance to a variety of problems 
in physical chemistry and biophysics [3]. In this study, we investigate the interac- 

tions between methane molecules and their charged derivatives in water-filled cylindrical 
nanopores. The interactions between these moieties are assessed in terms of the potential 
of mean force (PMF), which is the reversible work required to bring the solutes from infin- 
ity to r, the distance between the solutes. We are interested in such a geometry because, 
quantitative calculations of well defined systems might provide insights into more complex 
situations, such as the conformations of a polypeptide chain inside the cylindrical tunnel of 
the ribosome (4] and peptide transport through biological channels such as a-hemolysin jlj]. 
Our results are also applicable to other situations where the conformations of biopolymers 
are important, such as transport of nucleic acids through carbon nanotubes 6j and synthetic 
nanopores Q], and in nanopore sequencing of nucleic acids 8, 9]. Other situations of interest 
may include studies of adsorbed solutes in porous materials like zeolites. 

In a previous study we showed that hydrophobic interactions in spherical, nanometer- 
sized water droplets lOj differ significantly from those in bulk water. The PMF is inde- 
pendent of the size of the confining droplet and, in contrast to the bulk case, shows only 
a single minimum at contact; the solvent separated minimum (SSM) is completely absent. 
Interactions between oppositely charged solutes in spherical droplets show that the extent 
of solvation is determined by the charge density of the ion 10] . Here, we investigate the ef- 
fect of confinement in cylindrical nanopores on the interactions between methane molecules 
and their charged derivatives [llj. We show that, in a nanopore filled with water at the 
same density (p) as in the bulk, the SSM in the interaction free energy between methanes is 
absent. The contact minimum is metastable with respect to two other minima, at separa- 
tions corresponding to the diameter or the length of the pore. To understand the effect of 
confinement in a cylindrical volume on charged solutes, we also calculate the free energies 
of interaction between the ions M q+ and M g _. One of the methanes has a positive charge 
q + while the other carries a negative charge. The results illustrate how the balance between 
solute-solute and solute-solvent interactions depends on the charge density and the degree of 
confinement. We conclude by discussing the implications of our results for protein stability 
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in confinement. 



Methods: 

Models. We investigate the interactions between two methane molecules in cylindrical, 
water-filled nanopores using Metropolis Monte Carlo [3] simulations. We use the TIP3P 



model for water 



131 ]. and a unified atom representation for methane. The Lennard- Jones 
(LJ) interaction parameters involving the solutes [l4j are given in Tabled) Model ions M q+ 
and M _ are created by adding charges of magnitude q + = I c/ I = 0.4e or e to the methane 

n 

spheres, where e is the electronic charge [15]. 

Simulation Details. Simulations are performed at constant number of molecules N, 
total volume V and temperature T = 298 K, and thus sample the canonical ensemble. We 
assume that the effective volume available to the water molecules is 

V cS = V- N S V S (1) 

where N s = 2 is the number of solute molecules in the pore, and V s = (4/3)7r<r| f0 , is the 
excluded volume due to each solute. For a chosen number of water molecules, N w , and pore 
diameter D, the pore length L is calculated using 

V* = «(j\ L (2) 

and 

N w Tfl w . , 

Pw = 77 (3) 

VeS 

where p w = 997 kg/m 3 is the water density at 298 K and 1 atm pressure, and m w is the 
molecular weight of water. The pore dimensions and the corresponding N w are listed in 
Table [Til The potential energy due to the pore walls in the cylindrical polar coordinates £ 
and z is given by 

' i<D/2 
m) = < 12 ' and (4) 

10 12 kJ/mol £ > D/2 

,0 \z\< L/2 

U 2 (z) = { 11 " 1 (5) 

10 12 kJ/mol \z\ > L/2 

In order to simulate confinement effects we do not use periodic boundaries. Electrostatic 
and Lennard- Jones interactions are evaluated without a cutoff. 
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Monte Carlo trial moves consisted of a random translation for both solute and solvent 
molecules and a random rotation for solvent molecules only. In addition, we also used 



axis), of the methanes. The moves were accepted according to the Metropolis criterion |l2j . 
We used umbrella sampling to calculate the free energy of interaction between two methane 
molecules, —ksT log P(r), where fee is Boltzmann's constant and P(r) is the probability of 
finding the two solutes a distance r apart. Because of the quasi one-dimensional nature of 
the confinement when r ^> D, we do not subtract the free energy contribution — 2fc#Tlogr 
which arises from the increase in phase space proportional to r 2 in spherically symmetric 
systems. Therefore, these profiles cannot be compared directly to PMFs in bulk. The biasing 
potentials are harmonic, Ub = 0.5k(r — r ) 2 , where r is the center of each window and the 
spring constant k is chosen to be 5 kcal/(mol-A 2 ). Window centers are 0.1 nm apart in all 
calculations. The data is unbiased and the free energies are calculated using the weighted 
histogram analysis method (WHAM) [jij with code from Crouzy et al. fl7l |. All PMFs and 
free energies in this work are only determined to within an additive constant; the free energy 
at contact is arbitrarily chosen to be zero. 

Results and Discussion: 

Equilibrium structure of confined water. Fig. [Th shows the local solvent density 
divided by the average density p = N w /Veg for a pore with D = 0.8 nm, and number of 
water molecules N w varying from 10 to 128. For N w = 10 — 32, the reduced density is 
less than 1 for all values of £, implying that the excluded volume due to the solutes has 
been overestimated. For N w = 10 the solvent has no structure. As N w and L increase, two 
concentric water layers become progressively well defined. 

Fig. [lb shows the radial density profiles in nanopores with varying D and L and N w = 128. 
For the narrowest pores with D = 0.8 and 1.2 nm, two sharply defined concentric layers are 
seen. The volume available to the waters in the D = 0.8 pore is close to the volume accessible 
to water molecules in the periodically replicated channel with repulsive walls and diameter 



0.9 nm simulated by Lynden-Bell and Rasaiah [18j; correspondingly, solvent structure is 
similar in both cases. With increasing D, the layering is lost and the solvent becomes more 
bulk-like in the interior with substantial depletion at the surface (£ ~ -D/2). 

In Fig. [2^l we plot the probability distributions, P(/i 2 //i) for pores with D = 0.8 nm and 



collective translations and collective flips, z 



z (the z-axis being parallel to the por e 




4 



varying N w . The dipole moment of a single water molecule is fi, and \i z is its value along the 
z-axis. Fig. [2j3 shows P(fi z / aO for N w = 128 and varying values of D. In Fig. [2^, individual 
water dipoles are seen to be preferentially aligned with the pore axis, pointing either "up" 
or "down" . The extent of ordering increases with increasing pore length. Fig. [2b shows 
that, as the diameter increases from 0.8 nm to 2.0 nm, the strong ordering in water dipoles 
decreases and is altogether absent for the widest pore with D = 2 nm. However, there is no 
net dipole moment in the pores. The dashed lines show the probability distributions for the 
total dipole moment of all water molecules in a nanopore with D = 0.8 nm and N w = 128 
(Fig. Eh) and D = 2.0 nm and N w = 128 (Fig. [2b). The distribution in other nanopores is 
similar (data not shown). Thus, Figs. [T] and [2] show that solvent structure deviates the most 
from bulk behavior in long and narrow pores. 

Free energies of interaction between methane molecules. Contours of methane 
positional probabilities in nanopores, with N w = 128 and D = 0.8 nm in Fig. [3^ and for 
D = 2.0 nm in Fig. [3b show that the solutes are most likely to be found at the periphery 
of the pore where water hydrogen bonds are broken and water density is significantly lower 
(Fig. [lb). The effective interaction between these nonpolar solutes reflects the nature of the 
confining volume. Fig. H] shows free energies of interaction of two methane molecules in pores 
with D = 0.8 nm and L e [1-27, 2.59] nm. In all these pores, the solutes are most stable in 
contact with each other at r ~ 0.4 nm. The SSM that is seen in the methane-methane PMF 
in bulk water [l9,[2(3], is conspicuously absent. In pores with D = 0.8 nm, a distant minimum 
occurs at r m L. In all pores with L £ [1.27, 2.59] nm (Fig. HJ), the distant minimum is less 
favorable than the contact minimum by about 3 kJ/mol (~ 1.2k bT at T = 298 K). The 
height of the barrier between the two states increases from ~ 3.2 kJ/mol to ~ 8 kJ/mol 
with increasing pore length in this range, due to the greater free energy cost of disrupting 
the water column in the longer pores that are more ordered (Figs. [T]and[2]). In Fig. Owe 
plot free energy profiles (solid lines) in nanopores with D e [1.2, 2.0] nm and N w = 128. In 
all these pores, the contact minimum occurs at r = 0.4 nm. A distant minimum at r « 
max{D,L} (see Table HT1) . is marginally more stable than the contact minimum by about 
0.5 — 0.7 kJ/mol. There is another distant minimum at r w min{D, L}. In the case of the 
narrowest pores (D = 0.8 nm), the minimum corresponding to D almost merges with the 
one at contact and appears as a shoulder at 0.7 nm. With increasing D, the barrier between 
the contact minimum and the distant minimum at r ~ min{D, L} decreases, because the 
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greater water depletion at the pore surface (Fig. [Ib) lowers the free energy barriers and 
makes the states with intermediate r less unfavorable. 

Solute translational entropy due to surface solvation is a major contribution 
to the free energy. The free energy profiles in Fig. H] and can be qualitatively understood 
on the basis of the surface solvation of nonpolar solutes in confinement. At the hydrophobic 
pore boundary, the structure of water is disrupted due to the breakage of hydrogen bonds. 
This gives rise to a depletion zone of low solvent density at the boundary. At equilibrium, 
hydrophobic species such as methane will be localized at these extended cavities at the water 
surface. This result is further supported by quantum mechanical simulations which find that 
even an alanine zwitterion in an isolated water droplet 21[ is preferentially solvated at the 
surface, with the nonpolar side chain exposed to air. The entropic preference of nonpolar 



y replicated, cylindrical channels has also 
. The SSM seen in the PMF between two 



solutes for the surface of water-filled, periodical 
been reported in other computational studies jl8 ] 
methane-sized hydrophobic solutes in bulk water can be attributed to the presence of a 
single water molecule between the solutes, which is hydrogen bonded to other waters. In the 
solvent depleted zone near the walls of the pore where water hydrogen bonds are broken, 
this conformation is strongly unfavorable. Therefore, in hydrophobic confinement the SSM 
is absent. 

The contribution of the solute translational entropy to the interaction free energies can 
be calculated approximately, assuming that the hydrophobic species are hard spheres that 
are strictly confined to the surface of the nanopores. The translational entropy at a solute 
separation r, is proportional to the length of a contour defined by the intersection of a 
sphere of radius r and a capped cylinder with the dimensions of the pore, with the sphere 
center located on the cylinder boundary (see Appendix and Fig. [6]). The contribution of this 
purely geometric term to the methane-methane free energy, — TAS A , is plotted in black for 
three different pores in Fig. [5j As with the total free energy, — TAS A is only determined 
to within an additive constant, and is vertically aligned with the first distant minimum 
in the free energy at r « min{D,L}. Fig. [5] shows that —TAS A qualitatively captures 
the broad curvature of the free energy profiles, especially for the wider pores. Clearly, the 
solute translational entropy favors the first distant minimum in the free energy, while the 
solvent drives the solutes to either minimize or maximize their separation. We note that the 
solutes are not localized exactly at the pore surface (Fig. [3]), due to attractive van der Waals 
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interactions with the solvent. The first distant minimum therefore occurs at r < min{D, L}. 

Interaction between charged solutes: Extreme confinement induces ion-pair 
formation. In spherical water droplets [10], the tendency of solutes for surface or interior 
solvation depends on their charge density (£) and the curvature of the confining surfaces. 
For a given solute size and charge magnitude, anions have a greater preference for the 
surface than cations. The tendency for surface solvation in periodically replicated cylindrical 
channels also correlates with £. In simulations of monovalent cations and anions in a channel 
of diameter 0.3 nm jl8|. smaller ions were found close to the axis of the channel, while larger 
ones with the same charge were more likely to be located closer to the wall. 

The reduced preference of charged solutes for surface solvation, as compared to nonpolar 
solutes, is reflected in the free energies of interaction between model ions M q+ and M g _ of 
charge magnitude q + = \q~\ = 0.4e (Fig. [7]). In all the pores studied here, there is a contact 
minimum at 0.37 nm and a SSM at ~ 0.67 nm, with a difference of 2 — 4 kJ/mol between the 
two states. Relative to the contact state, the solvent separated state becomes more stable 
with increasing pore length at D = 0.8 nm (Fig. [7h), and also with increasing D at fixed 
N w = 128 (Fig. Eh). 

Fig. [8] shows the same free energies evaluated for solute charge magnitudes of e. Com- 
paring Figs. [T^l and [H^i, we see that increasing the solute charge in pores of D = 0.8 nm 
destabilizes the SSM. Similar findings hold good for the nanopore with D = 1.2 nm (Figs. [7b 
and [8b). However, for the wider pores with D = 1.6 and 2.0 nm, the opposite trend is 
observed; increasing solute charge makes the SSM more stable relative to the contact mini- 
mum. 

In the solvent separated state, the two solutes are stabilized in individual hydration 
shells with a single water molecule between them that is hydrogen bonded to other water 
molecules. From M^-water radial distribution functions (data not shown), we infer that the 
first hydration shell around the solvent separated solute pair with charge magnitude 0.4e, 
has a diameter of ~ 1.6 nm. For the solutes with a charge magnitude of e, the corresponding 
hydration shell has a diameter of ~ 1.3 nm, with a second hydration shell of diameter « 1.7 
nm. Thus, pores with D < 1.2 nm are too narrow to accomodate even a single hydration 
shell around the solvent separated solute pair with a charge magnitude of e. The direct 
solute-solute coulombic attraction is therefore stronger, and thus in these pores, ion-pair 
formation is promoted, especially when the charge increases from 0.4e, to e. The opposite 
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happens in pores with D > 1.6 nm in which the solute-solvent interaction energy is more 
favorable. Clearly, as D — ► oo one recovers the bulk behavior in which the ions are fully 
hydrated. 

Conclusions: 

The mutual interaction between solute molecules confined in a cylindrical pore is strongly 
influenced by the presence of the confining boundary. Solute-solute interactions differ drasti- 
cally from those in bulk solvent. Hydrophobic interactions are altered more than electrostatic 
interactions. Nonpolar solutes, such as the methane-methane pair, are driven to the pore 
surface to maximize solvent entropy. Consequently, the solvent separated state that is seen 
in bulk water is absent in nanopores. Methane-methane free energy profiles show one con- 
tact minimum and two distant minima, at separations corresponding to the diameter and 
the length of the pore. The broad curvature of the free energy profiles is shown to be largely 
determined by the translational entropy of the solutes at the pore surface. The solvent 
drives the hydrophobic solutes either to contact or to the largest possible separation in the 
pore. The intervening minimum arises due to the greater solute entropy at r « min{D, L}. 
The stability of the distant minima, relative to the contact minimum, depends on the pore 
dimensions. 

Model ions M q+ and M g _, created by adding charges of magnitude q + = \q~\ = 0.4e 
or e to the methanes, have a more favorable energy of interaction with the solvent. These 
charged solutes are solvated in the pore interior, away from the surface, and their interaction 
free energy shows a SSM. For a given charge density, this state becomes more stable relative 
to the contact state, as the volume of confinement increases. Increasing the charge from 0.4e 
to e increases the stability of the solvent separated state for the pores with D > 1.2 nm, 
which can accomodate a complete hydration shell around the ion pair. In narrower pores, 
increasing the ion charge destabilizes the SSM. 

The interactions between charged species in cylindrical confinement can be rationalized 
using an interplay of two length scales. In the presence of charges, an additional length scale 
/ b = z 2 e 2 /ek B T (z is the valence and e is the dielectric constant) plays an important role. The 
balance between Ib and D determines the stability of the SSM. To a first approximation, the 
stability of the second minimum (Figs. Hand [H]) is determined by Ib/D. As Ib/ 1 D decreases 
at a fixed charge magnitude, the second minimum becomes more pronounced. 
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The preference of aliphatic groups to be localized at the pore boundary and charged 
species to be solvated in the pore interior can be used to infer the stability changes of pep- 
tides confined to cylindrical pores. For a generic amphiphilic sequence it is likely that the 
reduction in the conformational entropy of the denatured states, with respect to the bulk, 
and the propensity of polar and charged residues to be fully solvated, should compensate 
for the tendency of hydrophobic species to be pinned at the interface. Thus, we predict 
that, confinement in water-filled cylindrical nanopores should enhance the stability of con- 



fined amphiphilic peptides. Preliminary molecular dynamics [22| and Langevin dynamics 
4| simulations support these arguments. The interactions between charged species M q+ and 
M q _ with q + = e suggest that ion pair formation, which is unlikely in the bulk, can be 
promoted in narrow cylindrical pores. The simulations (Figs. [7] and [S]) show that the in- 
teraction between charged solutes is governed by competing solute-solute and solute-solvent 



interactions, as well as the charge-density of the ions 10j. The importance of charge density 



in controlling the stability of folded RNA has been demonstrated previously 23]. 

Appendix: Estimate of entropy of pinning near the pore surface 

In this Appendix, we evaluate the free energy cost of localizing a spherical hydrophobic 
species near the surface of the pore. We idealize the hydrophobic solutes as point hard 
spheres that are strictly confined to the surface of the cylindrical pore. The free energy of 
interaction of these particles will be purely entropic, and can be calculated from the pore 
geometry. Consider an ideal particle located at the periphery of the pore, either along the 
curved length of the cylinder, or one of the flat end caps. The probability of finding the 
other particle at a distance r is proportional to the length of the contour defined by, the 
intersection of a cylindrical shell with the dimensions of the pore, and a spherical shell of 
radius r which is centered on the first particle (see Fig. [6]). Here, we derive the expression 
for this contour length. 

The equation for the body of a cylinder with its center at the origin and aligned along 
the z-axis is 

x 2 + y 2 = a 2 = (D/2) 2 , -L/2<z<L/2. (6) 

Here a, D and L are the radius, diameter and length of the cylinder respectively. The 
cylinder has caps at z — L/2 and z = —L/2. Eq. ([6]) can be parameterized as x = asin^ 
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and y = a cos 6. The problem is evaluated in three parts. 



I. SPHERE CENTER ON THE BODY OF THE CYLINDER, INTERSECTING 
THE BODY OF THE CYLINDER 

The equation of a sphere with its center at (a, 0, — c) on the body of the cylinder is 

(x - a) 2 + y 2 + {z + c) 2 = r 2 , -L/2 < c < L/2 . (7) 
Using the parameterized cylinder equations and solving for z in Eq. ((?]) yields z = — c ± 



y r 2 — 2a? + 2a 2 sin(6 l ). The length of the contour of intersection of the cylinder and the 
spheres centered on the rim of the cylinder at z = — c is 

h = Ana / y/(dx/d9) 2 + {dy/d9) 2 + (dz/d9) 2 dO, -L/2+c < ^r 2 -2a 2 + 2a 2 sin(0) < c+L/2 
Jo 

(8) 

The center of the sphere can be varied along the z-axis from -L/2 to L/2. Dividing by 2 to 
avoid double counting, the length of this part of the total contour is 

L/2 
-L/2 



Li(r) = 2na 



\ 1+ 2 f 2 Sm2 } T • /m d9 > - L / 2 + c ^ V / r 2 -2a 2 + 2a 2 sin(0) < c + L/2 1 dc 
Jo V r 2 — 2a 2 + 2a 2 sm{0) 



(9) 



II. SPHERE CENTER ON THE CYLINDER CAP, INTERSECTING THE BODY 
OF THE CYLINDER 

The equation of a sphere with its center on the cap of the cylinder at z — —L/2 and in 
the xz plane is 

(x - r'f + y 2 + (z + ^) 2 = r 2 , < r' < a . (10) 
Solving for z using the parameterized equations of the cylinder yields z = -| ± 



yV 2 — a 2 — r' 2 + 2ar'sin(#). The length of the contour is given by 

/•27T 

h = 4 / ^(dx/dO) 2 + {dy/dO) 2 + (dz/dO) 2 d6, < ^r 2 - a 2 - r' 2 + 2ar'sm(9) < L 
Jo 

(11) 
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The center of the sphere can be any where on the circle on the cap of the cylinder 
described by z = — L/2 and x 2 + y 2 = r' 2 . This gives a multiplicative factor of 2tit' . The 
center of the sphere on the cap can also be varied along the x-axis from to a. Therefore, 
the length of this part of the total contour is 




-5 r,2 ^ n (A) dg o < v / r 2_ a 2_ r ,2 + 2 ar>sm(9) < L \ dr' 

r l — a 1 — r l + 2ar sin (6*) 

(12) 

III. SPHERE CENTER ON ONE OF THE CAPS OF THE CYLINDER AND IN- 
TERSECTS THE SECOND CAP 

The equation of a sphere with its center on the z = —L/2 cap of the cylinder and in the 
xz plane is given by Eq. (110j) . Since the sphere should intersect the other cap of the cylinder 
present at z — L/2, substituting z = L/2 in Eq. ( flOi) . we get (x — r') 2 + y 2 = r 2 — L 2 , < 
r' < a. The equation for the cap of the cylinder at z = L/2 is 

x 2 + y 2 = a 2 , z = L/2. (13) 

The center of the sphere can be any where on the circle on the cap of the cylinder described 
by z = —L/2 and x 2 + y 2 = r' 2 . The center of the sphere on the cap can also be varied along 
the x-axis from to a. Therefore, the length of this part of the total contour is given by 

r a , | 7r ^Jr 2 - L 2 + r' < a, 

Ls(r) = / ^r'V^T 2 dr', 6=1 ( , (14) 

Jo cos" 1 [ r 2 + J( r i-L*) ) Vr 2 - L 2 + r' > a 

Finally, the total length of the contour of interest is 

l c (r)=L 1 (r)+L 2 (r)+L 3 (r), (15) 

where Li, L 2 anc 



L 3 are given by Eqs. (|T2l) and <HM respectively, and are solved using 
Mathematica |24j]. The probability of finding the two particles at a separation r is P(r) oc 
l c {r). The entropic contribution to the interaction free energy of the two hydrophobic solutes 
in this study, can then be approximated (to within an additive constant) as 

-TAS A (r)^-k B T\ogl c (r) , (16) 
11 



where ks is Boltzmann's constant and T is the temperature. 
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TABLE I: Lennard-Jones a parameters between the various sites. 





e a p [kJ/mol] 


cr Q/ 3 [nm] 


methane- methane 

oxygen-oxygen 
methane-oxygen b 


1.23 {€mm) 
0.64 (e 00 ) 
0.887 (ejuo) 


0.373 {(J M m) 
0.315 (a 00 ) 
0.344 (a MO ) 



Q The Lennard-Jones potential has the generic form V^j(r) = Ae a p (-^r-) 1 ^ — where a and (5 



identify the two species, r is the separation distance and a a p is the collision diprieter. 
fc The methane-oxygen parameters are obtained from the methane-methane 1J] and oxygen-oxygen 



3 



parameters by applying the Lorentz-Berthelot mixing rules. 



TABLE II: Nanopore dimensions for a given number of water molecules N w and number of solute 
molecules N m = 2. 



AT a 


D [nm] b 


L [nm] c 


10 


0.8 


1.27 


17 


0.8 


1.69 


26 


0.8 


2.23 


32 


0.8 


2.59 


128 


0.8 


8.31 


1.2 


3.69 


1.6 


2.08 


2.0 


1.33 



"Number of water molecules 
^Diameter of the cylinder 
c Length of the cylinder 
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Figure Captions: 

Fig. [TJ (color online) Radial density profiles for water in nanopores (a) with diameter 
D = 0.8 nm and number of water molecules N w varying from 10 to 128 (b) with N w = 128 
and D varying from 0.8 to 2.0 nm. The length of each nanopore is given in Table ITU 

Fig. [2J (color online) Probability distributions for [i x /fi, where /i is the dipole moment of 
a single water molecule and fi z is its component parallel to the pore axis, (a) with diameter 
D = 0.8 nm and number of water molecules N w varying from 10 to 128 (b) with N w = 128 
and D varying from 0.8 to 2.0 nm. The dashed lines show the probability distributions for 
the total dipole moment of all water molecules in a nanopore with (a) D = 0.8 nm and 
N w = 128 and (b) D = 2.0 nm and N w = 128. Water molecules are more aligned with the 
pore axis in long, narrow pores, but the net dipole moment is zero in all cases. 

Fig. EJ (color online) Contour plots of methane occupancies in nanopores with N w = 128 
and (a) D = 0.8 nm and (b) D = 2.0 nm. The occupancy probabilities of the nonpolar 
solutes are greatest at the periphery of the pore. The scale for probability densities is shown 
on the right. 

Fig. HJ (color online) Free energies of interaction of two methane molecules in nanopores 
with diameter D = 0.8 nm and number of water molecules N w varying from 10 to 32. All 
curves are translated vertically so that the reference state is at contact. The contact and 
one of the distant minima at r « L — a mm are clearly resolved. The other distant minimum 
appears as a shoulder at ~ 0.7 nm. 

Fig. (color online) Interaction free energies of two methane molecules in nanopores 
with N w = 128 and pore diameter D varying from 1.2 to 2.0 nm. The black lines show 
the contribution from the translational entropy of the hydrophobes, assuming that they are 
strictly confined to the surface. 

Fig. EJ (color online) Illustration of the intersection of a cylinder and a sphere with its 
center on the body of the cylinder. 

Fig. EJ (color online) Free energies of interaction of model ions M q+ and M g _ of charge 
magnitude 0.4e in nanopores (a) with diameter D = 0.8 nm and number of water molecules 
N w varying from 17 to 128 (b) with iV„, = 128 and D varying from 0.8 to 2.0 nm. 

Fig. [8j (color online) Free energies of interaction of model ions M q+ and M q _ of charge 
magnitude e in nanopores (a) with diameter D = 0.8 nm and number of water molecules 
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FIG. 2: 
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FIG. 6 
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